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Thermodynamic consistency of the Mean Spherical Approximation as well as the 
Self-Consistent Ornstein-Zernike Approximation (SCOZA) with the virial route to 
thermodynamics is analyzed in terms of renormalized 7-ordering. For continuum 
fluids this suggests the addition of a short-range contribution to the usual SCOZA 
direct correlation function, and the shift of the adjustable parameter from the po- 
tential term to this new term. The range of this contribution is fixed by imposing 
consistency with the virial route at the critical point. Comparison of the results of 
our theory for the hard-core Yukawa potential with simulation data show very good 
agreement for cases where the liquid-vapor transition is stable or not too far into the 
metastable region with respect to the solid state. In the latter case for extremely 
short-ranged interactions discrepancies arise. 
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I. INTRODUCTION 

By use of the Self- Consistent Ornstein-Zernike Approximation (SCOZA), very 
accurate results have been obtained for the equation of state for fluids^ and lattice 
systems (or the Ising spin system)^ 1 ^. Also in the critical region the results are 
very accurate, showing non-classical critical behavior-, but a form of generalized 
scaling is obtained in the critical region instead of full scaling 7 . 

Looking more closely into the data one finds that for the lattice gas, or nearest 
neighbor Ising model, results are somewhat better than those for the continuum 
fluid. The SCOZA critical temperature T c of the former is about 0.2 % away from 
the best estimates 8 , while this separation is about 0.6 % for the commonly consid- 
ered Yukawa fluid (z = 1.8) according to recent estimates based on Monte Carlo 
(MC) computations^. Here we will analyze this situation in view of renormal- 
ized 7-orderingi^iii, and we will propose a modified version of SCOZA. A slightly 
simplified variant of this scheme is implemented and solved numerically, yielding 
essentially perfect agreement with simulations for potential ranges relevant for the 
liquid state. For extremely short-ranged potentials where the gas-liquid transition 
is deeply buried inside the metastable region, however, agreement is less satisfac- 
tory, and the numerical procedure may even fail to converge to consistency with 
the virial pressure. 

As SCOZA builds upon the Mean Spherical Approximation (MSA), we can un- 
derstand that SCOZA is more accurate for lattice gases than for fluids by noting 
that this is true for the MSA, too. By further analysis it is seen that for fluids a 
significant contribution to the direct correlation function of short range is missing. 
When this contribution is included, it enlarges the amplitude of the correlation 
function and thus amplifies the internal energy due to correlations (which is neg- 
ative). This further increases the deviation from the corresponding mean field or 
Van der Waals type equation of state. Since correlations lower the critical temper- 
ature, the net result is a further lowering of it. For the Yukawa interaction typically 
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used as a model potential in simulations, we expect this additional lowering to be 
something like (0.5 ± 0.4) % as a crude estimate. An immediate problem here is 
the determination of the proper range of this short range piece of the direct corre- 
lation function. In this respect we can be guided by considering thermodynamic 
consistency with the virial theorem to estimate the range mentioned. 

In section |TT] we will investigate the relation between the MSA and renormal- 
ized 7-ordering for potentials of long range. Then in section II I II thermodynamic 
consistency with the virial theorem is considered and applied to SCOZA. A mod- 
ified theory is then formulated on this basis in section IIVI and applied to the 
one- Yukawa problem numerically in section |V] A short summary (section IVIj) 
concludes our contribution. 

II. RENORMALIZED 7-ORDERING AND MSA 

The MSA is the solution of the Ornstein-Zernike (OZ) equation, 

**> = t^wv 

with boundary conditions 

h(r) = — 1 for r < 1, 

(1) 

c(r) = —(3ip(r) for r > 1. 

Here h(r) is the pair correlation function, c(r) is the direct correlation function, 
and 

^,( r ) = -e-'tr-V/r (r > 1), 

is the perturbing interaction; the hard core diameter and the unit of energy are 
set to unity for simplicity; and a tilde denotes Fourier transformation. (It should 
be noted that ip(r) is not restricted to interactions of Yukawa form.) 



4 



For a lattice gas or an Ising spin system the boundary conditions become 

h(0) = -1, 

c(f) = —f3ip (r) for r 7^0. 

So in this case we have 



p + ?h{0) = p(l - p) = / - ^ d 3 fc (2) 



(2tt) 3 J 1 - pc(fc) 

for lattice cells of unit volume. (To keep the notation simple, we do not explicitly 
indicate that k and similar arguments are necessarily vectors in the lattice case, 
nor do we explicitly point out restriction to the Brillouin zone, i. e., —it < ki < tc 
with i = x,y, z.) Defining v(k) = —j3tp(k), we can now write 

c(k) = c + v(k) 

where cq is determined by Eq. (TSJ). For the reference system alone we then have 
$(k) = by which h(0) = h(k)) 



fi = p + p 2 h(k) = p(l -p) = (3) 

1 - pc 



More generally for ip(k) ^ we can define 

P 



c 



1 - pc ' 

and the core condition (EI) becomes 



fJL = p c + p\K 



(4) 



5 



K = —\ m ^ d"k. 



with 

(27r) 3 7 l-ti e v(k)' 
With n given by Eq. relation (j4j) can also be written as 

H = fji e {l- ^n"fJL e K), (5) 

where a prime denotes differentiation with respect to p. In this form we can 
recognize the relation obtained by the resummation or renormalization process to 
obtain the equation of state for fluids when the inverse range of the attractive 
interaction 7 is used as a perturbing parameter—^ (see Eqs. (26) and (36) of 



Ref. Illl ). In Eq. (jSJ), p c represents the renormalized hypervertex \x when graph 
expansion in terms of 7-ordering is considered^ 1 ^. The fi follows from the reference 
system, e. g., hard spheres, by use of the compressibility relation as 

/i = p(^)" = P + A(0), (6) 

where po and ho(r) are the pressure and the correlation function of the reference 
system. Thus /1 represents the Mayer graphs of the correlation function of the 
reference system. For small 7 they are regarded as <5-functions of relative distance 
r so that only the integral (jSJ) is needed as a leading order approximation. 

Now with renormalization (J3]) the equation of state (EOS) becomes more con- 
sistent, with better agreement between the energy and compressibility routes and 
essentially common critical point via the two routes^. But the remaining ther- 
modynamic inconsistency prevents well-defined isotherms near the critical point, 
which is also the situation of the MSA for the lattice gas or Ising model. 

In contrast to the MSA for lattice systems, the MSA for continuum fluids is 
not consistent with Eq. (jSJ) and therefore farther away from thermodynamic self- 
consistency. In fact, for interactions of very long range, %. e., small 7, the MSA for 
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continuum fluids also fulfills Eq. (j3J) but not Eq. (jSJ). This is due to the difference 
in reference systems. For a Percus-Yevick hard sphere reference system one has 



;i + 2 V y 

instead of expression ([3D, where 77 = |p for spheres of unit diameter. From this 
fi" = d 2 fi/dp 2 can be evaluated at the critical density, i. e., p c ~ 0.3 or r\ c w 0.15 
for interactions of realistic range. (In mean field, r] c = 0.129.) With 77 = 0.15 we 
find 

-H" = -0.53... » -0.5 (7) 

which is only one half of the "MSA value" that follows from Eq. (j4j). 

Compared with mean field the finite range of interaction lowers the critical 
temperature T c , and we have 

Tc « ^T C MF , (8) 
H 

where T C MF is the mean field value. (The spinodal curve and thus T c follow from 
the denominator of the integrand of integral (pO) being zero at k — 0.) In view of 
this it is clear that the MSA compressibility critical point is lowered too much so 
that the more accurate energy route gives a classical critical point well above the 
former. For interactions of shorter range where details of the core condition and 
hard core correlations become more important the above results will be modified, 
but the qualitative features of the MSA inconsistency will remain. 



III. SCOZA EVALUATIONS 



The SCOZA differs from the MSA in its imposition of thermodynamic consis- 
tency between the energy and compressibility routes to thermodynamics: A state- 
dependent free parameter is introduced into the pair structure and determined in 
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such a way that the consistency condition 

da ( 
d(3 = P 'dp 



da d 2 pu 

P^TJT' ( 9 ) 



u=-p I i[)(r)(h(r) + l)d 3 r (10) 



holds. Here 

7/ = 

2' 

is the configurational part of the internal energy per particle as computed from 
the energy route, and 

o=l-pc(0), (11) 

is the reduced inverse compressibility as obtained from the compressibility route. 

In earlier SCOZA evaluations the MSA form (TjQ) of the direct correlation func- 
tion has been used with (3 replaced by an effective value j3 e . The contribution to 
the internal energy per particle U\ from correlations is then given by integral (j3J) 

as 

- PUi = T^f^cK (12) 

with v(k) = —/3 e i/j(k). (Note that a "self-energy" term —p^{0) is included in (|T2"]) . 
but it will not contribute in ([9]) anyway.) In addition one has the mean field term 
Mo = p 2 ip(0)/2 such that u — Uq + U\. For continuum fluids p c should also depend 
upon k, i. e., p c = p c (k), which we disregard here to simplify the argument. Now 
at the critical point 1 — /i c /3 ec (— ip(0)) = is valid both for MSA and SCOZA that 
builds directly upon MSA via Eq. JT|. For the former (3 e = (3 so that (3 ec = /3 C , 
but one notes that both MSA and SCOZA have the same configurational internal 
energy at the critical point just as in the lattice case. 

Now with relation (jSJ) the internal energy u\ can no longer be equal to the MSA 
one at the critical point since the vertex function p c is different. According to the 
relation 

i - pM-^o)) = o 
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a change in /i c means a change in (3 ec while /j, c K in Eq. f|T2|) remains unchanged. 
The change in \i c thus has the consequence that the internal energy contribution 
from u\ at the critical point, u\ c oc l//3 ec = /x c , is modified, too. 

In SCOZA computations the shift in T c relative to the mean field value T C MF 
is determined by (JSJ) to leading order and not by the MSA /i c from (jlj) since the 
former is closer to thermodynamic consistency. At the same time the SCOZA 
energy U\ is dictated by the MSA // c . Thus there results an inaccuracy in the 
SCOZA u\ with a corresponding inaccuracy in the SCOZA T c . The reason is that 
to leading order the shift down in T c relative to T C MF is proportional to u\. In 
terms of Ui, Eq. (jSJ) is (/3 e = (3) 

fi = /i c (l + A), A = -//'/WP- ( 13 ) 

However, by solving SCOZA the tti is changed into the MSA value that follows 
from the fjt c = fif SA of Eq. ® instead of Eqs. (|5D or (fT3l). 

„MSA 
< SA =^m. 

He 

This new value modifies the equation (JT3l) into 

/ m msa\ 
/i = /i c f 1 + J . 

When inserted into ©, this gives an additional shift in T c compared with use of 
Eq. (USD, viz. , 
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Now with Eq. (JJJ) we approximate 

p^pf SA (l + p c K)^pf SA (l + 2A), 



where the last relation follows from comparison with Eq. (jSJ) and using the numer- 
ical value ((7j) for p". Inserted into Eq. ([H]) and expanded this means 

AT C « -A 2 T C MF (15) 

Evaluations for the Lennard- Jones (LJ) interaction gave^ 1 ^ 

u c T c 0.2545-4 

A = 1 - — = 1 = 1 = 0.064 

p T C MF 0.272 • 4 

so that 

AT C « -4.2 ■ 10" 3 T C MF » -4.5 • 10~ 3 T C . (16) 



IV. MODIFIED SCOZA DIRECT CORRELATION FUNCTION 

The lowering on T c given by Eq. ( fT5l) can be understood in terms of a graph 
expansion. There is a significant contribution of short range to h(r) and the direct 
correlation function c(r) that is not properly taken care of by previous SCOZA 
computations. To leading order it is connected to the //'-term of Eq. (jSJ). In terms 
of graphs, p represents the hypervertex where the graphs are those of the reference 
system correlation function plus S(f), i. e., p5(r) + p 2 h (r). On two p-vertices of 
these graphs a chain bond K with endpoints or vertices p, that can be replaced 
by p c to leading order, and symmetry factor ~ are added. This contribution is 
then added to the reference system correlation function or vertex p to obtain the 
renormalized vertex p^^-. 

For lattice gases, however, we do not need to consider this term explicitly as 
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the core condition (j3J) is then equivalent to (jSJ). To fulfill the former, a constant cq 
is added to the direct correlation function c(r) at r = 0. This can also be regarded 
as a short range piece of the perturbing potential such that v(k) in expression (jlj) 
is replaced by co + v(k). The consequence of this is that the chain bond integral (jlj) 
becomes zero and renormalization ([HJ) drops out as Co has already taken care of it. 

For continuum fluids, however, Eq. (jSJ) can not be replaced by a Cq(t) only 
inside hard cores if the core condition is to be kept; the short range part of it 
goes outside the core in contrast to the lattice gas system. Thus c(r) must have a 
contribution of short range outside the hard core (beyond that of the hard sphere 
reference system, chs(t)). This term will have both an amplitude and a range. 
In a SCOZA scheme this means two unknown parameters if now the remaining 
part of c(r) is just its MSA form outside the core, i. e., the new term replaces the 
earlier need for an effective temperature. 

Assuming the added short-range piece to the direct correlation to be of Yukawa 
form, too, the SCOZA direct correlation function for the one- Yukawa potential 
now becomes 

e -2i (r-1) g— Z2(r-1) 

c (r)=K x + K 2 + c HS (r), r > 1. (17) 

The first term can be treated as in MSA, %. e., its amplitude can be identified 
with the inverse temperature, K\ = (3, while its form coincides with that of the 
interaction, z\ = z\ this also means that c(r) asymptotically coincides with —/3ip(r) 
for large r as it should except at the critical point. K 2 and z 2 , on the other hand, 
are the two unknown parameters of the added short range piece to c(r). Obviously, 
conventional SCOZA is recovered when setting z 2 = Z\ = z, with K\ + K 2 playing 
the role of the effective inverse temperature f3 e . 

The use of an effective parameter j3 e is common to conventional SCOZA prob- 
lems. So the differences between various SCOZA problems lie in the actual equa- 
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tions used. These again depend upon whether continuum fluids or lattice gases 
are considered, and further they depend upon the pair interaction used. However, 
in the present case this conventional scheme is modified by which we have two un- 
known parameters present in Eq. (1171) . As only one parameter can be determined 
in "standard" SCOZA, an option is to keep the range fixed and to make evalu- 
ations for various values of z 2 to see how the estimated change in T c comes out. 
Furthermore, SCOZA allows the contact value of h(r) at the hard core surface to 
be obtained so that we can also check consistency with the virial theorem. Alter- 
natively, the range parameter z 2 can also be adjusted alongside K 2 as a function of 
(3 and p so that consistency between all three of the energy, virial, and compress- 
ibility routes is obtained throughout the temperature-density plane. In this more 
ambitious SCOZA evaluation first proposed by H0ye and Stell 14 both unknown pa- 
rameters, K 2 and z 2 , are supposed to be determined simultaneously. This scheme 
thus bears some resemblance to the version of the Generalized Mean Spherical 
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that did not, however, take into account cns( r ) 



Approximation studied in Ref. 
outside the core. 

Yet another possibility is to keep z 2 fixed for all densities and temperatures, 
choosing it in such a way that consistency with the virial theorem is achieved at the 
critical point. This has the advantage over full consistency of being much easier to 
implement while it can be expected that this simpler version already incorporates 
most of the effect of full consistency with the virial route, at least when using 
the location of the critical point and the phase boundary to gauge the relative 
accuracy and performance of various liquid state theories as we will do here. 
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V. NUMERICAL RESULTS 

In our computations we found it advantageous to adopt the Waisman recipe for 
the direct correlation function cns{ r ) of the hard sphere reference system, i. e., 

e -2ffS (r-l) 

CHs(r) = K H s , r > 1, 

r 

with density dependent parameters Khs and z hs that reproduce the Carnahan 
Starling equation of state 2 - 1 ^ 6 -. Just as in conventional SCOZA, the direct cor- 
relation function outside the hard core is then just a sum of Yukawa terms for 
which the solution of the OZ equation together with the core condition, h(r) = — 1 
for r < 1, is known semi-analytically by way of a Wiener- Hopf factorization and 
Eqs. ( TTU1) and (ITTj) can be evaluated efficiently^ 1 ^. 

The main change required vis-a-vis the conventional theory is the shift from 
Ki, which acquires an explicit temperature dependence, to K 2 as the adjustable 
parameter; also, it becomes necessary to implement the virial route expressions 26 
that do not normally feature in SCOZA computations but already appeared in the 
preliminary evaluation of our recent adaptation of SCOZA to molecules with soft 

19 

cores . 

For an interaction composed of a hard core of unit diameter with a perturbing 
tail ip(r), the virial theorem, 

/DO 
r 3 (h(r) + l)ip'(r)dr , 

relates the pressure P both to the contact value h(l+) of the pair correlation 
function and to an averaged inter-particle force away from contact. In the case 
of a direct correlation function of multi- Yukawa form outside the core, the same 
steps that lead to the energy route u and the compressibility route a, Eqs. ffTUj) 
and (1TT1) . also provide us with explicit expressions for the Laplace transform G(s) 



BP 2tt 

— = 1 + irP 
P 3 
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of r + r h(r)H from which the contact value can be worked out^. If, on the other 
hand, the interaction is a sum of Yukawa terms with inverse length scale parameters 
Zi, it is easily seen that the integral on the right hand side of the relation above 
reduces to a sum of linear combinations of G(s) and dG(s)/ds with s set to each of 
the z%. As both conditions are met in our case, the virial pressure can be evaluated 
numerically no less efficiently than the other thermodynamic quantities commonly 
entering SCOZA. 

Implementation of the theory outlined in section IIVI started from our re-imple- 
mentation of conventional SCOZA, which has been developed with a view to great 
conceptual simplicity, modularity, and flexibility and has turned out to be a con- 
venient testbed for modifications of both the theory and its numerical realization. 
Methodically it largely follows the work of Scholl-Paschinger— ; the main differences 
are the avoidance of general non-linear solving routines through the systematic use 
of linearizations (which also eliminates the need for cumbersome implicit differ- 
entiations and incidentally brings about a significant speed advantage) and an 
improved implementation of the artificial spinodal boundary condition (v. i.) that 
is more efficient computationally and does not restrict the estimated spinodal to 
the discretization grid, nor does it use the ill-founded criterion of minimal width 
of the spinodal. 

For the computations below, we solved the SCOZA PDE (JS]) by an iterated 
predictor-corrector method on density grids of N p = 10 3 density steps spanning 
the range < p < p max = 1 with the ideal gas (u = 0) and the MSA (K 2 = 0) as 
boundary conditions at p = and p = p max , respectively. The inverse temperature 
P was turned on from (3 = 0, where the hard sphere reference fluid provided the ini- 
tial conditions, down to the critical point with temperature steps decreasing upon 
approaching the critical point in approximate proportionality to min p a/ a(/5, p). 
Below the critical point the temperature steps increased again in such a way that 
the shift of the spinodal per integration step was limited to less than one density 
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step. Below the critical point there is a region where the analytical solution of 
the OZ equation becomes invalid^, and even a fully numerical procedure^ that 
does not rely on the analytical result eventually runs into problems of numerical 
convergence^. Either way, it proves necessary to remove part of the solution from 
the PDE's domain through the imposition of an artificial boundary condition at 
the spinodal. In our implementation we do so by estimating its location at the new 
temperature from the solution of the PDE at the previous temperature on nearby 
grid points and forcing the internal energy there to be compatible with a = 0. 
Pressure and chemical potential are obtained by thermodynamic integration at 
fixed density, starting from the hard core reference fluid at /3 = 0. This allows the 
binodal to be found by locating densities p v {0) and pi{0) of the coexisting vapour 
and liquid phases of equal pressure, P(p v ,f3) = P(pi,j3), and chemical potential, 
p(p v ,j3) = p(pi,P). The critical point is identified with the locus where the liquid 
and vapour branches of the spinodal and the binodal meet. 

This computation is then repeated for different values of z%. To illustrate the 
dependency of the critical parameters on z 2 , let us consider the potential with 
z = 1.8 that is known to exhibit thermodynamics roughly equivalent to the LJ 
interaction. As can be seen from Fig. [H the critical temperature obtained with 
z 2 > z i — z is indeed a bit lower than that for conventional SCOZA {z 2 = z\ = z), 
and typical shifts are comparable to the estimate ffTBj) obtained using results for 
the LJ interaction. 

As we solve the SCOZA PDE, we compare the pressure obtained from ther- 
modynamic integration to the virial route result. In general, for every value of z 2 
there will be lines along which consistency between SCOZA and the virial theo- 
rem occurs; this is illustrated, again for the case of z = 1.8, in Fig. [2J As can 
be seen, the geometry of these lines is rather complex, especially below the crit- 
ical point. We expect this distribution of z 2 values to resemble quite closely the 
function z 2 (P, p) that would be obtained with the more general version of the the- 
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ory imposing consistency with the virial route throughout the PDE's domain, cf. 
section ITVl 

In the present contribution, however, we consider a simpler version of SCOZA 
where z 2 is determined from consistency at the critical point only. Numerically the 
search for this optimal value of z 2 calls for repeated solutions of the SCOZA PDE 
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and might thus seem to be exceedingly burdensome. In contrast to, e. g., Ref. 
where this is indeed a concern, however, we can use the semi-analytic solution of 
the OZ equation. In combination with a very small number of Yukawa terms in 
Eq. ( TTT1) the integration of the PDE down to the critical temperature is sufficiently 
fast so that the search for z 2 does not constitute a problem. 

In tableUwe list the results both of the version of SCOZA outlined in section HVl 
(z 2 > Z\ — z) and for conventional SCOZA {z 2 — z\ = z) for various potentials, 
and we compare them to moderately recent simulations: Monte Carlo (MC) results 
obtained with N = 108 particles for z G {3.9, 7, 25, 100}^; Gibbs Ensemble MC 
(GEMC) data for z G {1.8, 3, 4} (N = 108, 9 ) and for z G {3.9, 7} (N = 216,% as 
well as the results of Grand Canonical MC (GCMC) simulations in combination 
with Finite Size Scaling (FSS) analysis for z = 1.82. The same data are also 
displayed in Fig. [31 It should be noted that the liquid-gas transition considered 
here is metastable with respect to the solid for rather short-ranged potentials, viz., 
for z > 6^; in particular, the highest z values for which simulations are available, 
viz., z G {25, 100} correspond to interactions for which the fluid-fluid transition is 
buried deep inside the metastable region^. 

First limiting ourselves to moderately short-ranged potentials [z < 10), a com- 
parison of conventional SCOZA (z 2 = Z\ = z) and of our modified SCOZA pro- 
posed here (z 2 ^ Z\ = z) shows two general trends: As expected, the critical 
temperature is shifted to lower values; for the LJ value of z — 1.8 this difference 
is about 0.7% and thus well compatible with the estimate ( fl6i) . Furthermore, the 
shift in critical temperature is also accompanied by a small decrease of the critical 
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density. Both effects become more pronounced as z becomes larger. The available 
simulations confirm that the shift in T c is certainly in the right direction, and 
the agreement is excellent considering the uncertainties of the simulation studies 
that generally cite non-overlapping confidence intervals for the same interaction. 
For some values of z the prediction of our modified SCOZA falls right into the 
confidence interval of some study, as is the case for z — 1.8 and the sophisticated 
analysis of Ref. ^ or for z = 4 and the GEMC result of Ref. O. For the nearby 
value of z = 3.9 deemed relevant to the fullerene Cgo, on the other hand, two differ- 
ent studies^ 1 ^ agree that our T c is too high, if by significantly different amounts, 
whereas our prediction is right between the results of the same two references for 
z — 7. For z = 3, our modified SCOZA is marginally compatible with the GEMC 
results. By way of contrast, conventional SCOZA (z 2 = z\ = z) yields too high 
a critical temperature in all cases. This situation is summarized graphically in 
Fig. [31 to judge by T c , the variant of SCOZA studied in the present contribution 
thus constitutes a genuine progress over its conventional form, and the available 
simulations are not sufficiently accurate to point out systematic deficiencies in 
the critical temperatures predicted. As for the critical density, both old and new 
SCOZA are well compatible with the results for z G {3,3.9,4}. Only for z = 1.8 
are we able to differentiate between the two versions of the theory: the modified 
SCOZA is in perfect agreement with the the extremely narrow confidence interval 
obtained by way of GCMC and FSS methods whereas conventional SCOZA lies 
somewhat outside this range, thus hinting at the superiority of our modifications 
for the critical density, too. 

As the potential becomes more and more short-ranged, however, discrepancies 
start to appear: For one, our modifications always bring about a shift of the 
critical density to lower values, but at z = 7 the GEMC p c is already higher 
than that of conventional SCOZA. Similarly, the critical temperature is always 
lowered, but the MC simulations for the two highest values of z listed in Tab. QJ 
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viz., z G {25, 100} yield critical temperatures that are far greater even than those 
of conventional SCOZA; furthermore, consistency with the virial theorem does not 
seem achievable at all for z = 100, presumably due to numerical difficulties with 
the evaluation of the virial pressure (cf. footnote I VII) . Clearly, more work is needed 
to clarify the situation at these extremely high values of z that are, however, far 
outside the range of stability of the liquid-vapour transition and thus of limited 
interest for a liquid state theory. 



VI. SUMMARY 



In the present contribution we have considered in some depth the question 
of thermodynamic consistency between MSA and SCOZA (that builds upon MSA 
but incorporates consistency between the energy and compressibility routes) on the 
one hand, and the virial route on the other hand. The basis for our investigations 
is formed by an analysis of the solution of the MSA closure in the presence of 
the core condition in the light of renormalized 7-ordering that highlights a crucial 
difference between the lattice and the continuum cases. For the latter it is seen 
that a short-range piece must be added to the direct correlation function, and it 
is this term that we propose to be adjusted instead of the usual potential term. 
This scheme has been implemented numerically in a simplified version where only 
the amplitude of this contribution to c(r) is varied throughout the p)-plane 
according to the usual SCOZA procedure implementing consistency between the 
energy and compressibility routes only whereas its range is fixed by imposing 
consistency with the virial route at the critical point only. Comparison with the 
results of various simulation studies shows that these modifications bring about 
significant improvements over conventional SCOZA; particularly encouraging is the 
perfect agreement with a sophisticated and extremely accurate study combining 
GCMC simulations with FSS analysis for z = 1.8 2 . 
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If all our evaluations have been restricted to the hard-core Yukawa class of 
interactions, the question arises naturally whether this approach can be applied to 
other types of potentials, too. In this respect it is clear that nothing prevents us 
from adding some Yukawa term to the MSA direct correlation function for, say, 
a Lennard- Jones potential and using the same strategies for evaluating its range 
parameter and solving the modified SCOZA PDE. At the same time, however, 
the form of the Yukawa potential is special in that it is most attractive right 
at the core whereas more realistic interactions reach the maximum of the depth 
of the potential only at greater distances. For these a contribution of Yukawa 
form to the direct correlation function hardly seems appropriate. So clearly more 
work is needed before the results presented here can be transferred to other cases. 
Nevertheless, we expect the strategy of using the virial theorem (that does not 
traditionally feature in SCOZA at all) as a means of gauging the fluid structure 
close to the repulsive core to remain the key for fixing the form of the adjustable 
part of the direct correlation function. This strategy, it should be noted, also 
featured prominently in our earlier work on SCOZA for molecules with soft cores 19 . 
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TABLE CAPTION 



TABLE 1: Comparison of the critical temperatures and densities for hard- 
core Yukawa fluids of various inverse range parameters z. SCOZA results were 
obtained as described in section |Vj with z 2 determined from consistency with the 
virial route as proposed in section HV] or set equal to z as in conventional SCOZA. 
The simulation results are taken from the literature. 



FIGURE CAPTIONS 



FIGURE 1: Variation of the critical temperature T c with the inverse range pa- 
rameter Z2 of the adjustable part of the direct correlation function for the potential 
with z = 1.8. The Z2 is varied from z 2 = z, corresponding to conventional SCOZA 
(T c = T c 5 ), up to z 2 = 600. 

FIGURE 2: Loci of consistency with the virial route for the hard-core Yukawa 
potential with z = 1.8 and various values of z 2 . Consistency is achieved at the 
critical point for z 2 = z 2 = 7.57145. The binodal obtained in this case is shown 
by the dotted line. This figure combines the results of several runs of our SCOZA 
implementation during each of which z 2 was kept fixed as opposed to a more 
general version of the theory where z 2 would be adjusted during a single run. 

FIGURE 3: Critical temperatures T c for hard-core Yukawa potentials with 
varying inverse length scale z as computed by conventional SCOZA (dashed line) 
and by various simulations relative to the prediction T* of our modified version of 
SCOZA. For the simulations, the crosses mark the value of T c , the vertical bars 
indicate the uncertainties where available, and the labels identify the source: "A" 
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- GCMC, FSS (Ref. a); "B" — GEMC 
Ref. 13); "D" — GEMC (N = 216, Ref. 



[N = 108, Ref. la); "C" — MC (N = 108, 



25|); cf. tabled 
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z 


T 


Pc 


method 


1 


2.51814 


0.279 


SCOZA, z 2 = z 




2.51388 


0.279 


SCOZA, z 2 = 5.95560 


1.8 


1.21869 


0.3145 


SCOZA, z 2 = z 




1.212(2) 


0.312(2) 


GCMC, FSS (Ref. 2) 




1.21013 


0.312 


SCOZA, z 2 = 7.57145 




1.177(5) 


0.313(13) 


GEMC (N = 108, Ref. 9) 


2 


1.08811 


0.323 


SCOZA, z 2 = z 




1.07877 


0.320 


SCOZA, z 2 = 8.03331 


3 


0.74025 


0.359 


SCOZA, z 2 = z 




0.72860 


0.352 


SCOZA, z 2 = 10.56167 




0.715(11) 


0.375(27) 


GEMC (N = 108, Ref. 9) 


3.9 


0.60200 


0.387 


SCOZA, z 2 = z 




0.58980 


0.373 


SCOZA, z 2 = 13.08509 




0.5714 


— 


MC (N = 108, Ref. 24) 




0.549(3) 


0.37(2) 


GEMC (N = 216, Ref. 25) 


4 


0.59108 


0.390 


SCOZA, z 2 = z 




0.57889 


0.375 


SCOZA, z 2 = 13.37749 




0.576(6) 


0.377(21) 


GEMC (N = 108, Ref. 9) 


5 


0.50885 


0.415 


SCOZA, z 2 = z 




0.49678 


0.392 


SCOZA, z 2 = 16.43066 


6 


0.45624 


0.437 


SCOZA, z 2 = z 




0.44390 


0.402 


SCOZA, z 2 = 19.75802 


7 


0.41881 


0.454 


SCOZA, z 2 = z 




0.411(2) 


0.50(2) 


GEMC (N = 216, Ref. 25) 




0.40530 


0.403 


SCOZA, z 2 = 23.44057 




0.4000 


— 


MC (N = 108, Ref. 24) 


8 


0.38977 


0.463 


SCOZA, z 2 = z 




0.37446 


0.397 


SCOZA, z 2 = 27.54452 


9 


0.36574 


0.465 


SCOZA, z 2 = z 




0.34842 


0.386 


SCOZA, z 2 = 32.09389 


10 


0.34497 


0.462 


SCOZA, z 2 = z 




0.32578 


0.372 


SCOZA, z 2 = 37.08103 


25 


0.2353 


— 


MC (N = 108, Ref. 24) 




0.18706 


0.307 


SCOZA, z 2 = z 




0.15829 


0.170 


SCOZA, z 2 = 165.80499 


100 


0.1538 




MC (N = 108, Ref. 24) 




0.06059 


0.07 


SCOZA, z 2 = z 
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0.999 



0.998 - 
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0.996 - 



0.995 - 



0.994 - 
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0.990 L 




Z2 
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